Raman transitions without adiabatic elimination: 
A simple and accurate treatment 



Rui HAN°, Hui Khoon NG a ' b and Berthold-Georg ENGLERT° ,C 

a Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore 
b DSO National Laboratories, 20 Science Park Drive, Singapore 118230, Singapore 
c Department of Physics, National University of Singapore, Singapore 117542, Singapore 

(Posted on the arXiv on 28 September 2012) 

Driven Raman processes — nearly resonant two-photon transitions through an intermediate state that is 
non-resonantly coupled and does not acquire a sizeable population — are commonly treated with a simplified 
description in which the intermediate state is removed by adiabatic elimination. While the adiabatic-elimination 
approximation is reliable when the detuning of the intermediate state is quite large, it cannot be trusted in 
other situations, and it does not allow one to estimate the population in the eliminated state. We introduce 
an alternative method that keeps all states in the description, without increasing the complexity by much. An 
integro-differential equation of Lippmann-Schwinger type generates a hierarchy of approximations, but very 
accurate results are already obtained in the lowest order. 
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1. Introduction 

Atomic and molecular systems can be coupled in various ways, so that the atomic states evolve 
and the populations of states change. Of all the electromagnetic multipole couplings, the electric 
dipole is the strongest. Thus optically allowed atomic transitions are often stimulated by an 
electric-dipole coupling between the atom and a well-controlled laser. 

Dipole-allowed transitions can be driven with optical lasers directly. It is, however, quite com- 
mon that the desired transition is dipole forbidden or the transition frequency is outside of 
the popular optical range. Then some intermediate state can assist in an indirect transition. 
Examples for such processes include three- level Raman transitions (1), multi- level Raman tran- 
sitions (0 and stimulated Raman adiabatic passage (STIRAP) Q |5j). They all have their 
useful applications, with three-level Raman transitions being perhaps the most widely used in a 
large variety of experiments. 

In order to reduce the computational complexity of dealing with large Hilbert spaces, one can 
decrease the dimensionality of the system by eliminating states that are not populated much 
or not coupled strongly. For a typical three-level Raman transition, the intermediate state is 
far off-resonantly coupled to the relevant initial and target states. This enables one to perform 
the so-called adiabatic elimination that gets rid of the less relevant state and yields a two-level 
effective Hamiltonian. Although the procedure of adiabatic elimination is well understood - 
see ([6j) and Q , for instance — it gives a reliable approximation only when the detuning of the 
intermediate state is much larger than the Rabi frequencies for the coupling to the other two 
states. 

In this paper, we introduce an alternative approach to the quantitative description of the 
three-level Raman transition. This new method does not rely on adiabatic elimination and gives 
a much more accurate solution while the computational complexity remains low. We set the stage 
in Section [2j where we briefly review driven three-level systems and state the notational conven- 
tions used throughout. Then, Section [3] deals with the usual adiabatic-elimination approach and 
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comments on its limitations and problems. Our new approach is explained in Section |4j First we 
present the general methodology, then we show how it gives the exact solution to the resonant 
two-photon transition problem, and finally we generate reliable approximations for the situation 
of a non-resonant two-photon transition. We close with a summary and outlook. 



2. The three-level system 

A Raman transition is a two-photon process that gives an effective coupling between two states 
|0) and |1) via a far-detuned auxiliary state |e); see Fig. [T] As mentioned above, Raman tran- 
sitions are often used when the transition between levels |0) and |1) is dipole forbidden or has 
an inconvenient frequency. The A-type configuration of Fig. [ija) applies to transitions between 
different ground states via an excited state; the cascade-type configuration of Fig. [jjb) can be 
used to achieve the transition between a ground state and a highly excited state, such as a Ryd- 
berg state. There is also a V-type transition with the level configuration similar to that of an 
inverted A-type transition. Although the level structures are different for these configurations, 
the underlying physics is essentially the same. 




Figure 1. Level scheme of a typical Raman transition, (a) shows the level structure of a A-type Raman transition and (b) 
shows the level structure of a cascade-type Raman transition. Qrj and fii denote the Rabi frequencies of the individual 
two-level transitions, A denotes the detuning of the laser from the transition frequency of the excited state and 5 is the 
detuning of the two-photon transition. The requirement is that the detuning A is much larger than the Rabi frequencies so 
that the excited state |e) is not significantly populated. 



We will, therefore, restrict ourselves to treating the A-type Raman transition in detail. With 
reference to |0), |1), and |e) in this order, the 3x3 matrix for the Hamiltonian of the system is 



H = H Atom + H AL , (1) 

where 

/0 \ 

H Atom = h\ Wl (2) 

\0 cjJ 
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is the part for the atom by itself with the convention that the state |0) has energy zero, and 

* / O e iWLO * \ 

Hal = J fi^' (3) 

1 ^e" iWL1 * / 

accounts for the atom-laser interaction. The laser frequencies are denoted by cjlo and cjli, 
respectively. The Rabi frequencies for the electric-dipole transitions are 

n = ^-(0\r-E LO \e) and ^ = ^(l|r • E u \e) , (4) 

where q e \ is the electron charge, -Exo and -Eli are the amplitudes of the electric fields of the laser 
beams and r is the position vector of the atom. Each Rabi frequency depends on the intensity 
and polarization of the corresponding laser as well as the dipole matrix element between the 
two coupled states, and can be complex. Typically, the detuning A is designed to be large, i.e., 
|A| 3> |Oo|, | so that the auxiliary state |e) does not get populated to avoid problems arising 
from uncontrolled spontaneous emission from |e) to other states; see Fig. [TJ Since the two-photon 
transition from |0) to |1) is nearly resonant, the overall detuning 6 of the two-photon transition 
is small compared with the average detuning A, |ff/A| "C 1. If the detunings between the lasers 
and the atomic frequencies are denoted by Aq = u e — wlo and Ai = lo c — u)\ — oju, we have 
A = (A + Ai)/2 and 6 = A - A 1 . 

In view of the time-dependent phase factors in Hal , it is expedient to switch to an interaction 
picture in which the Hamiltonian does not depend on time. This is achieved by identifying the 
"free" Hamiltonian Hq by a suitable splitting of the atomic Hamiltonian, 

±6 \ 

#0 = #Atom + M -iff \=h\ cut- ±6 ), (5) 
-A/ 




for which we obtain the interaction-picture Hamiltonian 



-6 fi, 



o 



Hi = e iHot/h (H - H )e- iHot/h = -\ 6 fii ] . (6) 



Q* 0* 2A 



The evolution of the three-level system is then studied with the aid of the Schrodinger equation 



ffij^l(t) = ftfl(i) with [ Cl (t) | -v ul "' ''*(/). (7) 



'co(t)' 

Cl(t) 
C e (t) 



where the components of ^j(t) are the interaction-picture probability amplitudes for |0), |1), 
and |e), related to the respective components of ^(t), the amplitudes in the Schrodinger picture, 
by simple time-dependent phase factors. As a consequence of this simple relation between the 
components of ^(t) and ^i(t), we can simply square co(t), ci(t), or c e (t) to obtain the probability 
amplitudes for the respective atomic levels, as exemplified by the probability for |0), 

10 0)^(t)| 2 = |(l0 0)^(t)| 2 = |c (t)| 2 , (8) 

where, of course, (l 0) is the three-component row for (0|. 

To understand the system analytically, we can solve for the eigensystem of this time- 
independent Hamiltonian directly. However, very often, this is neither the most efficient way 
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of getting the solution, nor the best method for obtaining a good physical insight into the sys- 
tem. Various approaches have been developed, of which the adiabatic elimination and, for 6 = 0, 
the dark-state method (0) are particularly useful and popular. 



3. Adiabatic elimination 
3.1. The methodology 

The standard textbook approach to the Raman-transition problem makes use of "adiabatic 
elimination" in accordance with the following line of reasoning. Since the excited state |e) is far- 
detuned by A, it will remain barely populated if it has no initial population. Thus, the change 

d 

of the population in this state can be taken as approximately zero, —c e (t) = 0, so that 

2i^c e (t) = n* c (t) + J2Jci(i) + 2Ac c (t) = 0, (9) 

in view of the Schrodinger equation Q. We can now express c e (t) as a linear combination of 
co(t) and ci(t), and so eliminate c e (t) from the equations of motion for co(i) and c\{t). This gives 
us an effective 2x2 Hamiltonian for the evolution of the two relevant states 

^=-U 5+ Jt for '4t m=Hcsm with m =(cM)- (io) 

V 2A " "r 2A / 

As an immediate benefit of applying adiabatic elimination on the intermediate auxiliary state, 
the effective Hamiltonian is a simple 2x2 matrix, for which the eigenvalues and the projectors 
to the eigenspaces are readily available. 
The eigenvalues of H e g are 

£ t = ~(|n | a + |n 1 | 8 )±^n H (n) 

with the positive frequency given by 

^ = (4A7^° |2 + |0l|a ) 2 + 4^° |2 - |fil ' 2 ) + 62 - (12) 

The projectors to the corresponding eigenspaces are |(1 ± a a ), where 

2H cS - E + - 

a ° = E + -E^ (13) 

is a Pauli-type matrix. If the evolution starts with all the population in the ground state |0), 
i.e., co(t = 0) = 1, the population in state |1) at a later time t is 

Ht)\ 2 = l -^^[l-cos(n R t)]. (14) 
R 

This tells us the physical significance of Or, = {E + — E_ )/h: It is the effective Rabi frequency 
of the transition between states |0) and |1) via this Raman process. 
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\ci(t)\ 



2| CIq I 2 | £2\ | 2 
l^oP + l^il 2 )' 



cos 



|^o| 2 + |^i| 2 
4A ' 



(15) 



where the effective Rabi frequency is Or = (|On| 2 + l^i| 2 ) /(4| A|) and the amplitude of the Rabi 
oscillation is less than unity, 



4|O | 2 |Oi| 2 
|^o| 2 + |^i| 2 ) 



< 1 



(16) 



unless |On| 2 = |Oi| 2 . In other words, when 5 = 0, we get complete population transfer between 
state |0) and state |1) only if the two lasers drive the respective transitions equally strongly. 

More generally, perfect population transfer from |0) to |1) is only possible if the two diagonal 
matrix elements of the effective Hamiltonian are identical, i.e., if there is no effective detuning 
after the adiabatic elimination. Then |0) and |1) are equal- weight superpositions of the eigen- 
states of H c s and temporal evolution turns one into the other. For given laser intensities, and 
thus given Rabi frequencies On and Oi, the experimenter can exploit the Zeeman or the Stark 
effect to adjust the overall detuning such that 



l^il 



I On 



4A 



(17) 



This makes the effective detuning vanish and ensures perfect population transfer. Indeed, for 



this value of <5, the right-hand side of (14) simplifies 



\ci(t)\ 



-[1 — cos (Or£)] with Or 



| On 1 1 0^ I 

2IAI 



(18) 



In experiments, specifically for two-photon population transfer from the ground state to a Ryd- 
berg state where the Rabi frequencies |On| an d can be an order of magnitude different in 



strengths, adjusting the detuning in accordance with (17) is important ©. 



3.2. Light shift 

When an atomic transition is driven by an electromagnetic radiation field with detuning A, the 
dressed atomic levels are shifted. This is the so-called light shift of the atomic levels, which is a 
second-order correction to the eigenenergies of the Hamiltonian (j KM . If an atomic transition is 
driven by Rabi frequency with detuning A, the light shift of the ground state of the transition 
is — 10| 2 / (4A) and the light shift of the dressed excited state is of the same amount but opposite 
in sign. The overall light shift is a direct summation of the light shifts arising from individual 
electromagnetic radiation fields, when the atomic level is addressed by multiple fields. 

In the present context, the light shift of state |0) is HSq = — fo|f2n| 2 /(4A n ) and the light shift 
of state |1) is HSi = — fo|Oi| 2 /(4Ai). Thus, the value of 8 that brings the two-photon transitions 
into resonance is determined by 
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Figure 2. Fidelit y be tween the two states at later time evolving with the exact Hamiltonian but for the two different 8 
values of fl7) and J20j. A = 400Mhz and = 40Mhz arc fixed and the different curves are for three different values of 
the ratio |s7fj/^i|- 



which is solved by 



2A(|^i| 2 - |fi | 2 ) 
8A 2 + |fi | 2 + l^il 2 



+ 



(20) 



where the ellipsis stands for terms of relative size |fi/A| 2 or smaller. The difference between the 
5 values obtained from (17) and (20) can be of the order of a percent (~ (|fi 2 + |^i| 2 )/(64A 2 )). 



But the approximate solutions provided by the adiabatic-elimination method do not depend 
much on this small fractional difference in 5. Figure [2] shows the fidelity between the states at 
later times using the two different values of 6 when the initial state is 1 0) . 

The effect of using these two different 5 values is small when the difference between |flo| and 
| fill is not too large. For the cases shown in Fig. § the error from using (|17|) is a small fraction 
of a percent. This explains why (17) can serve as a good guidance for experiments, although 



(20) is more accurate. 



3.3. Problems with adiabatic elimination 

In physics, one speaks of an "adiabatic process" if a relevant property evolves quite slowly, 
whereas other processes are fast — a clear separation of time scales is a defining element. In the 
context of adiabatic elimination, one invokes such a separation in the evolution of the components 
of v I / i(i), for which we have the Schrodinger equation in Q. One standard argument observes 
that, as a consequence of |A| S> |fio|> l^r |j the amplitude c e (t) will undergo many oscillations 
during a period in which co(t) and ci(t) do not change substantially. Then, if our interest is not 
in the short-time changes that would reveal the rapid oscillations of c e (i) but predominantly in 
the longer-time changes of co(t) and ci(t), the average change of c e (t) over these longer periods is 
expected to be quite small. In the spirit of this reasoning, we should then read ([9]) as a statement 
about coarse-grained values of the probability amplitudes. 

This is hardly a rigorous argument, and whether one regards it as convincing or not is largely 
a matter of taste. Clearly, though, a more solid argument would be welcome, and one has been 
provided in Q. Indeed, the reasoning in Q uses coarse graining in conjunction with a Markov 
approximation. 

There are other problems that one needs to keep in mind when employing adiabatic elimina- 
tion. We present them as four questions. 

(i) Which is the correct interaction picture to use? The theory argues that the change of the 
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(coarse-grained) population in the excited state is approximately zero and we arrived at 
([9]) by using the Schrodinger's equation of motion with the time-independent Hamiltonian 
H\ in ([6]). However, the choice of interaction picture is not unique and why do we apply 
the adiabatic elimination in this particular interaction picture instead of another? For 
example, by adding a constant term, the operator H\ + E describes the system equally 
well. But if we use an interaction picture with E 7^ 0, the resultant two-level effective 
Hamilton operator also changes as A gets replaced by A + E/h'm (10) with consequential 
changes in the evolution of ip(t). 

This ambiguity in the choice of interaction picture was studied by Brion, Pedersen, and 
M0lmer with a Green's function formalism (Tj5|). It would be premature, however, to claim 
that the case is closed. Further studies of the choice of interaction picture and a systematic 
way of performing the adiabatic elimination are presented in a companion paper Q . 

(ii) Is it possible to estimate the population in the eliminated state? This does not have a sim- 
ple answer because there are problems with the normalization of the wave function. Origi- 
nally, we have ^i(ty^i(t) = 1. Then, since H e g is hermitian, we also have ip(tyip(t) = 1. 
Combined with the initial condition (co,ci,c ) = (1,0,0) this implies c e (i) = for all t, 
which contradicts the basic approximation in ([9]). 

(iii) Is it possible to regard the adiabatic-elimination approximation as the first in a hierarchy 
of approximations? This is possible, indeed, as discussed in Q. It turns out that the next 
approximation in the hierarchy gives a substantial improvement and better quantitative 
estimates. In addition, it answers the previous question inasmuch as the next approxima- 
tion provides an estimate for the population in the auxiliary state. 

(iv) Is it possible to avoid the adiabatic elimination without increasing the complexity much 
beyond the convenient two-level description of (10)? Yes, this is possible, as we demon- 
strate in Section 01 



4. Without adiabatic elimination 

Alternative methods other than adiabatic elimination are also used to solve Raman-transition 
problems. The most direct way is, of course, to diagonalize the interaction-picture Hamiltonian 
of ([6]), and this can be done by hand since the dimensionality of the system is small. But the 
expressions for the eigenvalues and eigencolumns of Hi are quite involved and not transparent. 
An exception is the two-photon transition with vanishing overall detuning, 5 = 0, when one can 
identify a dark state and use it to reduce the three-level system to an effective two-level system, 
which can then be solved exactly rather simply. The dark states are particularly useful in the 
context of adiabatic population transfer and electromagnetically induced transparency. When 
the overall detuning is nonzero, however, there is no dark state. Methods of perturbation theory 
can then be used to find corrections for a small detuning 5, but the complexity grows quickly 
when high accuracy is required. 

In this section, we provide a new way of solving the Raman-transition problem. Just like the 
dark-state method, it gives a compact exact solution for 5 = 0, and it can solve the 5^0 case 
with very high precision and rather little extra effort. 



4.1. General methodology 

It will be expedient to use a different interaction picture as the one of ^ and ([6]) that we used 
for the adiabatic elimination in Section 3.1. Instead, we choose 



'A + 5 
H = £ I 2wi + A - 5 
2uj e - 



(21) 
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and the Hamiltonian in this interaction picture is 



A-5 


o si 





-A + 6 Sl! 


o* 

"0 


n; a 



where the latter way of writing emphasizes the 3 = 2 + 1 split into two relevant states and one 
auxiliary state; 03 is the standard third Pauli 2x2 matrix, and Si and Sl^ are the two-component 
column of Rabi frequencies, 

( 23 ) 



and its adjoint row. This change of interaction picture is equivalent to shifting the energy levels 
of the interaction-picture Hamiltonian in ^ by —HA/2, an example of the freedom of choice 



discussed in question (i) in Section 3.3 
We note that the square of Hi, 

H\ = {hM) 2 = h 2 (Mjj + e), (24) 
is the sum of a "big" block-diagonal part, 

Mo =4~V A 2 + f^J' ( 25 ) 

and a "small" off-diagonal part 

_ 5( a 3 Sl\ 

Matrix M denotes a square root of the matrix (Hi/h) 2 and, since for each eigenvalue we have a 
choice of sign, there are many Ms that are equally good. Any one can be used as a replacement of 
Hi/h in even functions of Hj; for example, we could choose M > by convention. In particular, 
then, the unitary evolution matrix in this interaction picture can be written as 

U(t) = e~ iH ^ h = cos(Mt) - '^M±± Hl . (27) 

Except for the common physical approximations that enter the modeling of the atom-laser system 
by a driven three- level system described by the Hamiltonian of Q-Q, this is an exact 3x3 
matrix representing the evolution operator. 

Instead of diagonalizing the interaction-picture Hamiltonian ( J22| ) , we can determine the eigen- 
values and eigencolumns of M 2 , whose "big plus small" structure, together with the block- 
diagonal form of Mq , facilitates approximations. We will see the advantage thereof shortly. Let 



us note that we can position the factor H\ in the second term of the right-hand-side of (27) 
equally well to the left of sin(Mt)/M, or break up sin(Mt)/M and sandwich H\ between even 
powers of M. Since [Hi, M 2 ] = 0, such a change in the order of the matrices makes no difference 



in (27), but slightly different expressions are obtained when approximations are introduced for 



the trigonometric functions of M. 
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4.2. Vanishing overall detuning (d — 0) — exact solution 

When 5 = 0, we have e = and 



M = Mi 



o 



VA 2 + Mt 

VA 2 + fitf] 



(28) 



where the signs are chosen such that M — > H\/h in the — > limit. The resulting evolution 
operator (27) reads 



U(t) = cos(M t) 



i sin(M t) 
h Mo 



(29) 



Owing to the block-diagonal structure of Mq, the original 3x3 problem has been converted 
into an equivalent 2x2 problem without introducing any approximation. Clearly, the technical 
difficulty has been significantly reduced! 
The column Q is an eigencolumn of A 2 + f2f^ with eigenvalue A 2 + Qtf2, so that 



M 2 = i(A 2 + Ot )^oj + Ia 2 ^ 











(30) 



is the spectral decomposition of Mq. This gives 

(OQt 
mi 



+ cos(At/2) 



not 



o 



mo 





(31) 



and likewise for sin(Mot) in ( |29| ). Hence, the exact evolution of the system can be written out 
analytically. In particular, the population in the excited state |e) is non-zero. With the atom 
initially in the ground state |0), it is 



(t)\ 



;ooi)c/(t) 



a 2 + ntn 



sinf\/A 2 + OtJ2t/2 



(32) 



This exact expression shows that the population in the excited state oscillates with (angular) 
frequency 7j(A 2 and the oscillation amplitude can be non-negligible if |fio| 2 i s a sizeable 



fraction of A 2 + = A 2 + 



Oil 



We note that the second projection matrix in (30) and (31) projects on the dark state, whose 
bra has the row (— Q\ ^Iq 0) (8). An atom prepared in this dark state stays in it, and there 
is no probability of finding the atom in the excited state at any time. The atom is essentially 



decoupled from the driving lasers under these circumstances. In this sense, one could regard (29) 



with (30) as the evolution matrix in the dark-state formalism but this is, in fact, not the case. 



In the dark-state approach, one diagonalizes Hi, reduced to a two-dimensional problem after 
putting the dark state aside, which amounts to choosing one particular square root of Mq from 
the continuous family of square roots that the degenerate eigenvalue makes available, namely 
the square root whose eigenvalues and eigencolumns are those of H\/h. No such unique Mq is 
needed in (29), nor is there any benefit from enforcing a unique square root of Mq by imposing 
additional criteria. Although in the 5 = case, this equivalence can be established between the 
new approach and the dark-state approach, the new approach offers more flexibility and it also 



has a clear advantage when dealing with the 5^0 case, as we shall see in Section 4.3 
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Figure 3. Population distribution for a single-atom Raman transition in time when 5 = 0. The solid curves show the exact 
solution, and the dashed curves the adiabatic-elimination approximation. The initial state is |0). The red curves are for the 
ground-state population |co(t)| 2 , the blue curves for |ci(t)| 2 , and the orange curves report |c e (t)| 2 , the population in the 
excited state. The detuning is A = 400 MHz for all plots; the top left plot is for |f2 | = l^ll = A/10 = 40 MHz; the top 
right plot is for |fl | = A/10 = 40 MHz and \Qi\ = A/16 = 25 MHz; the bottom plot is for |f2 | = l^ll = A/4 = 100 MHz. 



For three sets of parameter values, Fig. [3] shows the population |c e (t)| 2 of (32) as well as the 
populations |co(t)| 2 and |ci(i)| 2 of the relevant states |0) and |1). The populations of |0) and |1) 
are oscillating with a slow frequency that goes with a large amplitude and a fast frequency that 
goes with a small amplitude. The population transfer is predominantly controlled by the slow 
frequency which, therefore, is the effective Rabi frequency of the system. Since M5 2 has only two 
distinct eigenvalues, half the sum of their square roots gives the fast frequency and the difference 
gives the slow frequency. Thus, we have 



for the effective Rabi frequency. The expansion in powers of |fi/A| 1 permits a comparison 
of this exact expression with the approximation obtained by adiabatic elimination, the 5 = 



version of (12): The adiabatic-elimination approximation gives the leading term, but does not 
reproduce any of the higher-order terms. 

In summary, we find that the solution from the adiabatic elimination is indeed the zeroth-order 
approximation of the exact result in the expansion of £t*£l/ A 2 . For example, in the zeroth order 
of fi^O/A 2 , the excited state population vanishes as the oscillation amplitude is proportional 
to f2tf2/A 2 ~ 0. We conclude that adiabatic elimination yields a reliable approximation only 
when |Oo|i|^i| The difference between the solutions obtained by our method and the 

adiabatic-elimination approximation is demonstrated clearly in Fig. [3| 

The top left plot in Fig. [3] shows the populations when | Oo | = |^i| = A/10 is small; in this 
parameter regime, we can already see the deviation of the adiabatic elimination from the exact 
result but the deviation is not significant. In the top-right plot, the population of state |1) only 
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reaches a maximum of about 80% as |0q| 7^ |^i| ~ A/10 introduces an effective detuning for the 
two-photon transition; the adiabatic elimination gives a good approximation since we are still 
in the regime of \£Iq\, |^i| *C |A|. When the magnitudes of |flo| and are not much smaller 
than |A|, the population of the excited state is no longer negligible, but the complete population 
transfer between states |0) and |1) is still achievable; see the bottom-left plot in Fig. [3j The 
adiabatic elimination method fails in this case of stronger coupling between the relevant states 
|0), |1) and the excited state |e). 



4.3. Non-zero overall detuning (5 y£ 0) 



As discussed in Section 3.3 one needs to adjust the detuning 5 to compensate for the light shifts 
and achieve complete population transfer from |0) to |1). Besides this, there can also be other 
experimental reasons for choosing a particular 5 value. Thus, the situation of 5 7^ is of practical 
interest, and so we will regard the overall detuning 5 as a free parameter that is small compared 
with the average detuning A. Typically, the ratios 5/ A and O'JI/ A 2 are of the same small order, 
a few percent perhaps. 



The splitting of Hf in (24)— (26) has the diagonal blocks in Mq, including the <5-dependent 



contributions, whereas e contains the off-diagonal parts linear in 5. We shall treat e as a small 
quantity, without, however, regarding 5 as an expansion parameter. Rather, the full <5-dependence 
of Mq is taken into account. Then, the eigenvalues of M 2 agree with the eigenvalues of Mq to 
first order in the perturbation of e, and the leading correction will be of order e 2 . 
The eigenvalues of M5 2 in the subspace spanned by |0) and |1) are 



M | = l( A 2 + j2 ) + \tftt±\J (^Slf + %5A^<j& + {45A) 2 , (34) 
4 8 8 * 

which we get from diagonalizing the upper 2x2 diagonal block (A + 5a^) 2 + . The unnor- 
malized eigencolumns are 

/[V-(A + fo 3 ) 2 ]f!\ , ,_, 



V 



for M = /4 , (35) 



unless Q, is an eigencolumn of 03, which is a case of no interest. The eigenvalue | (A 2 + QJQj in 
the lxl block for |e) does not depend on 5. With its eigenvalues and eigencolumns at hand, all 
functions of Mq are readily evaluated. 
In passing, we observe that the difference between \xh and p?_ is smallest, as a function of 5, 



when (17) holds, i.e., 4<5A = -fiV 3 a Then 



/4 " M 2 - = \ V (^f " = ^ , (36) 

which is nonzero in all situations of interest. 

For a systematic inclusion of correction of orders e, e 2 , e 3 , . . . , we do not use the perturbation 
theory for an approximation of the eigenvalues and eigencolumns of M 2 for use in (27). Rather, 
we generate approximations for the evolution matrix U(t) = exp(—iHit/h) itself with the aid of 
an equation of Lippmann-Schwinger type, 

un-uPm-fv ^'w-vK uv), ( 3 7) 
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where 



U^\t) = cos(M t) 



i sin(M t) 
H W 



(38) 



differs from (29) by the inclusion of the 5 dependent terms in Mq and H\. One way of verifying 



that (37) is correct, is by checking that both sides have the same Laplace transform. Indeed, 



they do: 



1 



1 



dte s U(t) - s + - g2 + M 2 hg2 + 



s 2 + M 2 s + iHi/h 



(39) 



is an identity that follows from (24) 



As mentioned at the end of Section 



4.1 



the multiplication order of Hi and even powers of M 2 
is irrelevant in (27) as they commute. In (38), however, the order does matter as [Hi, Mq] ^ 
when 5 ^ 0. In addition to (37) with (38) where Hi is on the right, we have, therefore, also an 
"on the left" version, 



V^-tv sin(Mot ° 



o 



M 



(40) 



with 



(41) 



Half their sum gives a symmetrized version, 



u(t) = uP(t) 



dt' 



sin(M (t -t')) 
Mo 



with 



rA s )u\ (\/f+\ 1 sin(M t) i sin(M t) 

U Q (t)=cos(Mot)--^—Hi--Hi^— 



(43) 



and there are many more variants that one could explore. 



Each of the integral equations (37), (40), and (42) provides a hierarchy of approximations by 



an iteration that commences with the respective zeroth-order approximation. This is, of course, 
the procedure by which one generates the Born series from the Lippmann-Schwinger equation. 



For the "on the right" equation (37), the A;th-order approximation is 



U m {t) = U w {t) 



(R) 



dt 



, sin(M (t-tO) (R) 



M 



(44) 



and analogous expressions apply to the "on the left" version and the symmetrized variant. Note 



that Uj. (t) is not half the sum of Up(t) and Uf'(t) for k ^ 0; that arithmetic mean could 
also be taken as a valid fcth-order approximation. For such a scheme to be useful in practice, the 
zeroth-order approximation should be quite good to begin with, the first-order approximation 
should be sufficient for many purposes, and it should not be necessary to go beyond the second 
order. 

The various fcth-order approximations differ from each other, but they are all accurate up to 
fcth-order in e/Mg ~ |5S7|/A 2 . In the common case where |^/A| ~ 10 _1 and \S/ A\ ~ 10~ 2 , we 



r(R) 



(L), 
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Figure 4. Plots of populations obtained from different zeroth-order solutions of the Lippmann-Schwinger equations. Blue 
curves give the exact results from numerical simulation using the Hamiltonian Hj\ green curves show solutions from the 
symmetric approximation E/g (t); red curves show solutions from t/g R '(i); and orange curves show solutions from t/g L '(t). 
The parameters are A = 400 MHz, |fl | = A/2, = 3A/10 and 8 = -Stf <rz£l / {4A) = (\fti\ 2 - |£7 | 2 )/(4A) = -16MHz. 
The effective Rabi frequency is = 27.8 MHz, about 7% of A. Initially, we have co(0) = 1 and ci(0) = c e (0) = 0. The 

curves starting at 1 show the approximations for |crj(t)| 2 ; the curves that start at and rise to 1 are for |ci(t)| 2 ; and the 
curves that start at and never exceed small values are for | c c (t) | 2 . 



have |e/Mg| ~ 10~ 3 . Remember that the 5 dependence in Hi and the e dependence (which also 
depends on 5) are treated separately. Although e goes to zero when 5 vanishes, \e/Mq 2 \ is one 
order of magnitude smaller than |<5/A|. 

A technical point is the following. The approximate evolution matrices Uk(t) are not unitary, 
rather Uk(tyUk(t) deviates from the unit matrix by an amount of order e fc+1 . One can cope with 
this in various ways (1 1 Perhaps the simplest is to ensure proper normalization by including a 
time-dependent factor that depends on the initial set of probability amplitudes, thereby arriving 
at an effectively unitary matrix Uk(t) that is suitable for the given initial column ^(O), 

g>(t)*,(0)= , U J- mm (45) 

In other words, we apply Uk(t) to ^i(O) and normalize the resulting column to unit length. This 
procedure worked fine for all examples that we studied. 

As remarked above, the fcth-order approximations uj^ (i), uj^\t), and lljf\t) differ slightly 
and might not describe the system equally well. We discuss a few examples of the state popu- 
lations as functions of time under the different approximations, and analyze their performance. 
We compare the results to the exact numerical answers. The value A = 400 MHz is taken for 
the average detuning, as it is of typical order for real experiments, and the values of the Rabi 
coupling strengths fio and f2i can vary in a range of fractions of A. The overall detuning 5 of the 
two-photon transition can be controlled to within 1 MHz accuracy in laboratory experiments. 
For the purpose of this analysis, then, we take the liberty of setting 5 to any value we like. 

Figure [4] gives an example that demonstrates the quality of the zeroth-order approximations 
given by f7 (R) (i), U^\t), and U^(t). Since uP(t) is an average of U^ R \t) and U^\t), we 

expect the population curves for Uq (t) to lie between the other two curves and this can be 
seen quite clearly in Fig. |4j Moreover, the solutions of the population in the initial state |0) do 
not depend much on which of the three zeroth-order approximations is used, and all of them 
are very close to the exact numerical solution. In this example, we use 6 = (|^i| 2 — |Qq| 2 )/(4A), 



the value of (17), which gives an effective resonant two-photon transition, and we can see that 
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Figure 5. Comparison of the probabilit ies obtained from the zeroth-, first-, and second-order approximations of the 
Lippmann-Schwinger equations (37] and |4p| . In plot (a), the blue, green, and red curves show solutions from U^\t), 
[/£ '(t), and U2 [t), respectively; and in plot (b), the blue, green, and red curves show solutions from from C/g L '(t), 
(7-f L ^(t), and U^(t), respectively. The parameters and the initial state are the same as in Fig.^ 

full population transfer from |0) to |1) can be achieved. In comparison, full population transfer 
cannot be achieved for other values of 5. 

In Fig. El the value of \5\ is about a few percent of |A|. The plot shows that, in this parameter 
regime, Uq \i) approximates the evolution of the probabilities for finding |0) and |1) quite well, 
and it certainly works best among the three different zeroth-order approximations shown here. 
The approximation for the excited state population also works well when t is short, i.e., during 
the first few fast oscillation periods, but the deviation grows quickly with time. 

The accuracy is better for higher-order approximations. We compare the approximations of 
zeroth, first, and second order for the three different Lippmann-Schwinger equations in Figs.[5ja), 
[j^b), andj^j FigureJ^a) shows that the deviation of the zeroth-order approximation U^\t) from 
the exact numerical solution is large when about half of the effective Rabi cycle is completed, 
i.e., around t = 7t/T2r or At ~ 45, and the deviation is smaller around a full Rabi cycle. The 
main deviation is in the size of the small-amplitude oscillations with short period, whereas 
the Rabi oscillation with longer period is reproduced equally well by Uj^ (t) with k = 0, 1,2. 
The first-order approximation corrects part of the error in the zeroth-order approximation, and 
the second-order approximation improves matters further and gets the probabilities very close 
to their exact values. The same observations can be made about the corresponding "on the 
left" approximations in Fig. [s|b). How about the symmetric version uj^\t) whose zeroth-order 
approximation already works quite well? 

Figureplshows that when the two-photon transition is resonant, the difference between Ugit), 

— (s) — — (S) J I 

U{ (t), and U% (t) is difficult to detect. All three lowest-order approximations of (42) describe 



the evolution of the system well during the first Rabi cycle. The zeroth-order approximation 
works surprisingly well, at times it gives a better result than the higher-order approximations 
(see the bottom plot of the two blow-ups). The improvement offered by the higher-order approx- 
imations can be observed near the middle of the Rabi cycle (see the top plot of the blow-ups) , 
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Figure 6. Comparison of the zeroth-, first-, and second-order approximations of the symmetric Lippmann-Schwinger l |42[ i. 
The parameter values, the initial state, and the color coding are the same as in Fig. [5] 

and this improvement is more substantial when e gets larger. For a practical application to ex- 
periments that aim at complete population transfer from |0) to |1), the approximation provided 

~ (s) 

by Uq (t) is accurate enough to determine the parameter values reliably. 



4.4. Discussion 



To summarize our approach, we split (Hi/h) 2 , the square of the interaction Hamiltonian, into 
two parts: the dominant part Mq and a small correction e. With this splitting and any one of 
the Lippmann-Schwinger equations (37), (40), or (42), successive approximations to the unitary 
evolution matrix U(t) = exp(—iHit/h) can be obtained iteratively. If we use the approximations 



given by uf^\t) and Ul u> (t), we only need to do one iteration and stop at the first-order solution 

[k = 1) for a very good approximation. If we use the symmetric version f/£ (i), we do not even 

need to go beyond the zeroth-order approximation since Uq (t) is already very close to the exact 
evolution for typical experimental parameters. Thus, we have 



f(L) 



U(t) ~ U^ S \t) = cos(M f) 



i sin(M t) i sin(M t) fr( s )(+\ 



(46) 



where U, 

r(S) 



(S) 




(t) differs from U Q 



(S) 



(t) by the time-dependent factor of (45) that ensures *&i(t) 



Uq (t) x I / i(0) is properly normalized for the given initial column of probability amplitudes. 



According to (46), the oscillation of the state populations of |0) and |1) are governed by 



the operator Mq. The effective Rabi oscillation frequency, to a very good approximation, only 
depends on the eigenvalues of the first diagonal block of Mq. Applying the same argument 
as in Section 4.2 for the case of 5 = —Q^a^Q/AA in (17), the effective Rabi frequency is (we take 
M± > 0) 



2V( A2 + ^ + 2 ( 



inol + inil) 2 



ij(A2 + *») + i(|no| 



(47) 



With this particular choice of 5, the Rabi oscillation amplitude could reach unity; however, this 
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Exact 

uf\t) 

Adiabatic 
elimination 



Figure 7. Improvement on the effective two- level Hamiltonian compared with adiabatic elimination. The black curves 
give the exact numerical solution; the green curves are for the symmetric zeroth-order approximation of U(t); the blue 
dashed curves are for the adiabatic-elimination approximation; and the red curves result from taking Mo as the effective 
Hamiltonian. The parameters are the same as in Fig. [4] 

might not be the real maximum that the population in state |1) can reach, because on top of 
this slow effective Rabi oscillation, the population also oscillates with a fast frequency. This fast 
oscillation goes roughly with the frequency /z+ + fj,- ~ |A|. Nevertheless, up to linear order in 
8/ A, we find that the oscillation amplitude goes to unity, regardless whether we choose 5 in 
accordance with (17) or with (20). This shows why the evolution of the system is essentially the 



same for both 5 values. 

Moreover, since the effective Rabi oscillation of the two relevant states depends on the eigen- 
system of Mq only, the effective Hamiltonian between states |0) and |1) is approximately given 
by the 2x2 upper diagonal block of HMq, if we exclude the excited state |e) from the evolution 
directly. That is 



h 



(A + cuts) 2 + not . 



(48) 



The minus sign is chosen with the same reasoning as in (28). From this effective Hamiltonian 



we can find the Rabi frequency directly and one of the special cases was already given in (47) 
the oscillation amplitude for an arbitrary 5 is 



(!^m) 2 + 8£Afit<7 3 fi + (45A) 2 ' 



(49) 



and again we have P = 1 when 5 = — f^<73f2/4A. Figure [7] shows that this effective Hamiltonian 
is much more accurate than that of the adiabatic-elimination approximation, inasmuch as the 
evolution by e~ lMot gives a very close envelope of the population oscillation. 



5. Summary and Outlook 

After reviewing the standard adiabatic-elimination approximation, which reduces the theoreti- 
cal description of a three-level Raman transition to an effective two-level problem, and identi- 
fying some of the shortcomings of this approach, we introduced an alternative approximation 
method. Similarly to adiabatic elimination, there is an essential two-level component in the new 
method without, however, eliminating the third auxiliary level. This makes the new method 
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easy to use, inasmuch as one only needs to diagonalize a 2 x 2 matrix. Integro-differential equa- 
tions of Lippmann-Schwinger type are the powerful tools that enable us to generate successive 
approximations. A particular one with high symmetry performs so well that the lowest-order 
approximation is all one needs for a highly reliable determination of experimental parameters. 

For the sake of simplicity in presentation, the discussion was here limited to two-photon 
transitions. The method can also be applied to more complicated situations ([77]) . about which 
we will report on another occasion. 
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